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Abstract 

The grand potential of a classical Coulomb system has universal finite-size corrections 
similar to the ones which occur in the free energy of a simple critical system : the massless 
Gaussian field. Here, the Coulomb system is assumed to be confined by walls made of an 
ideal conductor material ; this choice corresponds to simple (Dirichlet) boundary conditions 
for the Gaussian field. For a d-dimensional [d > 2) Coulomb system confined in a slab of 
thickness W, the grand potential (in units of ksT) per unit area has the universal term 
T{d/2)({d) /2'^'n:^/'^W'^~^ . For a two-dimensional Coulomb system confined in a disk of 
radius R, the grand potential (in units of ksT) has the universal term (1/6) InR. These 
results, of general validity, are checked on two-dimensional solvable models. 
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1. INTRODUCTION 

A classical Coulomb system is a system of charged particles, interacting through the 
Coulomb potential, plus perhaps some short-range interaction ; there may also be a con- 
tinuous charged background. We are interested in equilibrium properties, and classical 
(i.e. non-quantum) statistical mechanics is used. A Coulomb system may have phase 
transitions (at least in the two-dimensional case) ; here we assume the system to be in a 
conducting phase. 

In some geometries, it has already been shown that the free energy (or the grand 
potential) of a classical Coulomb system exhibits universal finite-size corrections, in close 
analogy with what happens in critical systems^^"^^ . In Coulomb systems, although the 
screening effect makes the charge correlations short-ranged, that same screening effect 
makes the electric potential and field correlations long-ranged^^^ ; the critical-like behavior 
of the free energy is related to the existence of these long-ranged correlations. The electric 
potential is in some sense the analog of the order parameter of a critical system. 

For defining a Coulomb system with boundaries, one has to state the boundary con- 
ditions. In a previous paper*^^-', we considered the case when the boundary is a plain hard 
wall which confines the charges. Such an assumption does not generate any simple bound- 
ary conditions for the electric potential : for a given charge distribution, the potential 
"leaks out" , taking in general non-zero values outside the Coulomb system and that spoils 
the analogy with the order parameter of a critical system. 

The purpose of the present paper is to study classical Coulomb systems with bound- 
aries, in the simpler case of Dirichlet boundary conditions for the electric potential : the 
walls, impenetrable to the particles of the Coulomb system, are assumed to be made of 
an ideal conductor material on which the electric potential vanishes. In general, one must 
also assume that there is some wall-particle short-range repulsion which prevents a particle 
from collapsing onto its electric image. 
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Two geometries will be considered. 

(a) Slab. In d dimensions (d > 2), a Coulomb system is confined between two parallel 
ideal conductor plates, separated by a distance W ; the Coulomb system extends to infinity 
in the d — 1 other directions. Let uj be the grand potential per unit area. We shall show 
that uj (times the inverse temperature (3) has the large-14^ expansion 

/3a; = + — ^ + ... (1.1) 

The first two terms represent, respectively, the bulk and the surface contributions ; the 
coefficients A and B are non-universal. However, the last term of (1.1) is universal, with 
a coefficient C{d) depending only on the dimension d : 

where T and C are the gamma function and the Riemann zeta function. In particular, 
C(2) = 7r/24. 

Before we derived the expansion (1.1), we guessed it by analogy with a similar expan- 
sion which is valid for a simple critical system : the massless Gaussian field theory. We 
define the partition function of that theory as the functional integral on a field 0(r) 



Zg = y"l?0(r)exp ct>{v){-A)^{v)dv 



(1.3) 



where /i2 = 27r, (U3 = An, and more generally, for d > 2, = {d — 2)2n'^/'^ /T{d/2). Some 
ultraviolet regularization is needed for avoiding divergences. With Dirichlet boundary 
conditions for 0(r) on two parallel plates, it was shown^^) that the free energy associated 
to Zq, in the case d = 2, has a universal term —7r/2AW. More generally, we shall show 
that the expansion of the free energy for the Gaussian field theory is of the form (1.1), 
except for a change of sign of the universal coefficient C (d) . 
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(b) Disk. A two-dimensional Coulomb system is confined in a disk of radius R. The 
boundary circle is assumed to be an ideal conductor. For that case, we shall show that the 
grand potential Q has the large- i? expansion 

pn = AR^ + BR+-\nR+ (1.4) 
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where again the coefficients A and B of the bulk and perimeter contributions are non- 
universal, while the term (1/6) In is a universal finite-size correction. Here too, we first 
guessed (1.4) by analogy with a similar expansion^^'^^ which holds^ for the Gaussian field 
theory, except for a change of sign of the universal term. 

In two dimensions, there are exactly solvable models of Coulomb systems : the two- 
component plasma and the one-component plasma, at some special temperature. The 
general expansions (1.1) and (1.4) can be explicitly checked on these models. 

The critical system which is related to Coulomb systems is the massless Gaussian 
field theory ; this relation is discussed in Section 2. Section 3 is about the slab geometry. 
Section 4 about the disk geometry. The Conclusion suggests some further research. 

2. COULOMB SYSTEMS AND GAUSSIAN FIELD THEORY 

There are many signs that there is some connexion between Coulomb systems and the 
Gaussian field theory. The Coulomb potential is the inverse of the properly normalized 
Laplacian — (l///(i)A which appears in (1.3). If (f) is interpreted as the electric potential, 
the Hamiltonian in (1.3) 



is the familiar expression for the Coulomb energy in terms of the electric field — v</'(r). 
Beyond some microscopic cutoff distance, the correlations of the electric potential in a 

^ More general domain shapes were considered in refs. 5 and 6. 




(2.1) 



Coulomb system are Gaussian^^^ and of the same form as the correlations associated with 
(1.3). 

However the partition function of a Coulomb system is not identical to (1.3), in par- 
ticular because (1.3) involves a functional integral on the potential rather than the familiar 
integral on the particle positions. A Gaussian transformation was made in ref. 1 for relat- 
ing the Coulomb system to (1.3), but, although this has been useful for the slab geometry, 
we have not been able to deal with the disk geometry by this method. 

We suggest here another approach, based on the heuristic assumption that the uni- 
versal features of the grand partition function of a conducting Coulomb system are 
correctly accounted for by the functional integral expression 



fpexp 



£ 
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j dr dr' p(r) ^(r, r') p(r') 



(2.2) 



where p(r) is a charge density and G(r, r') the Coulomb interaction solution of 



AG'(r,r') = -//d 5{v - v') 



(2.3) 



with the appropriate (here Dirichlet) boundary conditions. In (2.2), the non- universal par- 
ticle structure of the Coulomb system is already disregarded, and some cutoff prescription 
has to be made. Performing now a change of function from the charge density p(r) to the 
electric potential 0(r), which are related to each other by 

A0(r) = -Hd p(r) 



we obtain 



Dp_ 
D(f) 



Zg 



(2.4) 



where Zq is the Gaussian partition function (1.3) and Dp/Dcp the Jacobian of the trans- 
formation from p to (p. Since 



Zg 



det ( — -A 



6 



and 




one obtains, from (2.4), Ec = I/Zq and 



In Sc = — In Zg 



(2.5) 



Thus, In Ec and In Zq have the same universal term, except for its sign. The change of 
sign comes from the Jacobian Dp/Dcf) in (2.4), i.e. the replacement of the charge density 
by the electric potential, as already explained in a different language in ref. 2. 



We consider a slab in d dimensions. Let us write the (/-dimensional position vector as 
r = (x, y), where y is the dth component and x stands for the d—1 other components. The 
Coulomb fluid fills the slab between y = and y = W, where there are ideal conductor 
walls. 

In this Section, we derive the large- expansion (1.1) in general, in two different ways, 
and we check its validity on two two-dimensional solvable models. We make a comparison 
with the case of plain hard walls. 

3.1. Gaussian field theory route 

For the free energy of the two-dimensional Gaussian field theory in a strip, with 
Dirichlet boundary conditions, the expansion (1.1) with C = — 7r/24 was derived in ref. 3. 
The relation with Coulomb systems, discussed in Section 2, gives for Coulomb systems the 
expansion (1.1) with C = 7r/24. 

These considerations can be generalized to d dimensions, for instance by adapting a 
method used in refs. 3 and 4, as follows. By considering the functional integral (1.3) as 
a path integral on an imaginary time interval W, one can express (1.3) in terms of the 



3. SLAB 



Hamiltonian of a quantum field theory in cZ — 1 space dimensions. This Hamiltonian is a 
sum of terms 

where k labels all possible wave-vectors in a space of dimension d — 1 and 11^ is the 
canonical momentum conjugate to 0k- Because of the Dirichlet boundary conditions = 
at "times" and W, (1.3) is expressed in terms of (0, 0) matrix elements : 

ZG = TT(0k = O|e-^^H0k = O> = nf • 

Therefore 

f d'^-^k ( (3k 



For substracting off the ultraviolet divergences, it is convenient to consider first 

df 1 r d'^-'^k 

(3—— = - / —k ctnh Wk 

dW 2 J (27r)rf-i 

and remove the W = 00 value, which gives 

1 r d'^-^k 



df df 
P—!- = (3- 



dW dW 



W=oo 



1 f d'^-'k , 



dW 



^ (d-i)r(f)c(rf) 



Thus, (3f has the universal finite-size correction —C{d)jW'^~^ for the Gaussian field 
theory, and (from Section 2) the correction is C{d)/W'^~^ for a Coulomb system, with 
C{d) given by (1.2). 

3.2. Screening sum rule route 

For the slab geometry, a more direct derivation of (1.1) is possible. It avoids any 
explicit reference to the Gaussian field theory, and therefore does not rely on the heuristic 
connexion of Section 2. Instead, this derivation uses a sum rule which expresses the perfect 
screening property of a conductor. 
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Let us first consider a cZ-dimensional Coulomb system which fills the half-space y > 0, 
with an ideal conductor wall at y — 0. Let Ey{0) be the microscopic electric field at 
some point, say the origin 0, on the boundary and p(r) the microscopic charge density at 
r = (x, y) (y > 0). The sum rule is 

p J dr y (^Ey{0)p{r)y = -1 (3.1) 

where < • • • >^ means a truncated equilibrium statistical average : < AB >^=< AB > 
— < A >< B >. One can derive*^^-' (3-1) by considering the response of the system 
to an external infinitesimal point dipole p oriented along the y-axis and located at the 
origin (on the Coulomb system side). This dipole adds to the Hamiltonian an interaction 
term —pEy{0), and by linear response theory the average charge density at r changes by 
5p{r) = f3p (^Ey{Q) p(r)^ . Now, we assume that the Coulomb system has good screening 
properties : the induced charge density (5p(r) is localized near the origin and has a dipole 
moment which cancels p : 

j dry 8p{v) = -p (3.2) 

(since the dipole moment is defined with the origin chosen on the wall, there are no 
contributions of the surface charges on the ideal conductor wall). The sum rule (3.1) 
follows. 

Let us now add a second ideal conductor wall at y = W . The good screening prop- 
erties imply that, if W is large, the local structure near the origin is unchanged, up to 
exponentially small effects. Therefore, the charge density 5p{r) is unchanged, its dipole 
moment (3.2) is unchanged (the integration range on y in (3.2) can be kept as (0, oo) since 
5p{y) is localized near y = 0), and the sum rule (3.1) remains valid for the two- wall system. 

This sum rule (3.1) can now be used for computing a universal correction to the 
pressure on the wall at y = 0. In free d-dimensional space, the Coulomb interaction, 
solution of (2.3), is 

G'o(r,r') = -ln(|r-r'|/a) , d^2 
Go(r,r') = l/|r- , d>2 
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(a is some irrelevant length scale). With the present boundary conditions, the solution of 
(2.3) becomes 

CO 

^(r, r') ^Go{r- r') - Go{r* - r') + [Go{r + n2Wu - r') 

n=l 

-Go{r* + n2Wu - r') + Go{r - n2Wu - r') - Go{r* - n2Wu - r')] (3.3) 

where r* = (x, —y) is the image of r = (x, y) with respect to the wall at y = and where 
u is the unit vector along the y axis. Using 



EyiO) = - I dr- ^ ' ' 



(0) = -/, 



dy' 

we can write 

T f dG{r,r') 



P(r) 

r'=0 



EyiOfY = -J dr 



dy' 



r'=0 



Since ^£^y(0)p(r)^ is short-ranged, we can expand dG/dy' in powers oiW at fixed r, 
with the result 



+0 (1/iy'^+i) (3.4) 

where e{d — 2) = d — 2iid>2 and e{d — 2) = 1 if d = 2. Let us describe the short-range 
interaction between the wall at y = and the particles by some potential F(y), and let 
n(j/) be the particle number density. The pressure is 



f°° dV 1 / ^ \ 



where the last term is the electrostatic pressure [if the short-range interaction is an infinite 
barrier at y = yo, —dV/dy must be replaced by f3~^5{y — yo)]. Since the screening effect 
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makes the density profiles near the wall y = independent of W (up to exponentially 
small corrections), the related average field 

oo 



is also independent of W . Thus, the pressure (3.5), where {Ey{QYj can be rewritten as 
(EyiQfY + (^y(O))^' depends on W only through (3.4). Using the sum rule (3.1) in 
(3.4) gives for the pressure p 

(d-i)r(|)c(rf) / 1 \ 

(3p = MW = oo + ^ ' /^^^^ ' + O —— . 3.6 

Since a;, the grand potential per unit area, has the derivative 

du) 

-P 



dW 

the grand-potential expansion is indeed given by (1.1) and (1.2). 
3.3. The two-dimensional two-component plasma 

This is a two-dimensional system of particles of charges ibg. The model is exactly 
solvable*^^'^°^ at a special temperature f3~^ such that f3q^ = 2. For f3q^ > 2, a point- 
particle model is unstable against the collapse of pairs of oppositely charged particles ; we 
work with almost-point particles, introducing some short-distance cutoff in the interaction. 
In presence of ideal conductor walls, for fiq^ > 1 the model would also be unstable against 
the collapse of a particle onto its electric image, unless some short-distance cutoff is also 
introduced in the wall-particle attraction. 

In two dimensions, it is convenient to represent a position r = (x, y) by the complex 
coordinate z — x + iy, and the Coulomb interaction (3.3) can be explicitly summed into 

sinh k{z — z') 



G{r, r') = - In 



sinh k(z — z') 



(3.7) 
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where k = 7r/2W (this result can also be obtained by a conformal transformation from the 
half- plane) . The calculation for the present strip geometry is very similar to the one which 
has already been done''^^-* for the half-plane case with only one ideal conductor wall. We 
start with a lattice model to avoid the collapse of positive and negative particles. There 
are two interwoven lattices and L_ with iV_|_ and N- sites. The complex coordinates 
of the sites of (resp. L_) are {ui}i<i<jv+ (resp. {fi}i<i<jv_ )• Positive particles can 
be on the sites of L+ and negative ones on those of L_. Each site contains at most one 
particle. We work in the grand canonical ensemble and denote the fugacity by (. The 
interaction between two particles is given by (3.7), and furthermore each particle has a 
self-interaction^ — (q'^/2) In |aA;/ sinh A;(2; — ^)|. Writing 



2 

and using exp(2A;«) and exp(2/cf) instead of u and v in the Cauchy double alternant 
formula, we can follow the same steps as in ref.ll and we obtain for the grand partition 
function S = det(l -|- K) where K is the matrix 



K 



/ f iakC \ f iakC \ \ 

y sinh A; (iti— it J ) J 1<»<JV+ ysiTihk{ui—Vj) J i<i<N_^_ 
l<j<JV_|_ 1<J<^V- 

/ —iakC, \ I —iakC, \ 

1 \syah.k{vi—Uj) J i<i<N_ ysinh k{vi—Vj) J i<i<N_ 

\ 1<J<-'V+ l<j<JV_ / 



We now consider the continuum limit where the lattice spacing goes to zero. We define 
the rescaled fugacity m = 27ra(/S where S is the area of a lattice cell. The eigenfunctions 
('0, x) of m~^K and its eigenvalues 1/A are defined by the two coupled integral equations 

/ =^^(^') + 7\X{^'^) = V'(r) (3.8a) 

4:W Jd smhk{z — z') smh A;(2; — 2;') 

-jk / ^'^'^^17^7^^^^^') + -^:Tr- — r^xi^'^) = x(r) (3.86) 

4:W Jd smhA;(2; — z') smhA;(2; — z') 



^ In equation (2.3c) of ref. 11, exponents 1/2 are missing for the first two factors on 
the r.h.s. 
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where the domain of integration D is the strip. In terms of these eigenvalues, 



A 

By using the equahty 

d 1 TT 



lnS = 5]ln(l + ^) . (3.9) 



5{r-r'), r,r'eD 



dz sinh k{z — z') k 

we can transform the integral equations (3.7) into the differential equations 

d d 
2^V' = ^I'X and 2— zx = 
oz oz 

In terms of the Dirac operator <t.V) these equations have the familiar^-'^^) form 

(<T.v)* = A* (3.10) 

where here is the spinor {ij^^ix}- These equations can be combined into the Laplacian 
eigenvalue problem 

AV' = A^V (3.11) 
with boundary conditions given by the integral equations (3.7): 

'4>{x, 0) = — x(a;, 0) and i^ix, W) = x{x, W) 

(incidentally, with such boundary conditions, —A turns out to be non-hermitian). We 
look for a solution of the form 

i;{x,y) = e'^'^'iAe'^y + Be-'^y) with I = iy/v"^ + A^ (3.12) 

The boundary conditions give two linear homogeneous equations for the coefficients A and 
B. The determinant of this system must be zero for a non- vanishing solution to exist. 
This gives the relation between v and A: 

sinh(W^V^;^ + A2) 

cosh{wVv^ + A2) - A \ = (3.13) 
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Let us define the entire fonction 

1 



cosa^Wy + z-^) — z — 



cosh f 

The solutions of (3.13) are the zeros of / and we have /(O) = 1, so 

n 

Ae/-i(0) 

Therefore, for each value of f , a partial summation of (3.9) can be made, giving for the 
grand potential u per unit length 

1 /■+°° ^ , 1 /•+~ 



1 /■+°° TO 1 /"+°° 

/?a; = / c^t'ln TT (1 + — ) = / dv\af{-m) 



Ae/-i(o) 

Now, in the limit W ^ oo ^ 



ln/(-m) ~ Wi^/v'^ + m2 - |i;|) + ln(l + ) - ln(l + g-^l'^l^) + O (e""^^) 

The first term gives the bulk contribution to the grand potential, the second the surface 
term and the third the universal finite-size correction. The bulk and surface terms would 
diverge if we did not take into account the particle-particle and particle-wall short-range 
repulsion. We can do so, for a repulsion of range cr, by introducing a cutoff Vmax = '^/cr. 
The final result for the grand potential per unit length can be written as 

(3u = -(3p,W + 2/37e + + 0(e--^) (3.14) 

where 

w? f 2 \ , ^ 

/?P6 = — In + 1 (3.15) 

ZTT \ ma J 



and 



m ( 2 7r\ 
/57c = -— In— + 1-- (3.16) 

l-K \ ma 2 



The bulk pressure is essentially the same as the one found for an infinite system ^^^^^ up 
to a slightly different cutoff prescription. The surface tension 7c for ideal conductor walls 
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is given by (3.16). And, finally, we find the finite-size correction C{2)/W = 7r/24VF, as 
expected. 

In the argument of Section 3.2, it was assumed that the local structure near one wall 
is unchanged, up to exponentially small corrections, by the presence of another wall at 
a large distance W. This assumption can be explicitly verified on the present model, by 
a calculation of the one-body and many-body densities, by the same method as in some 
previous work*^^°'^^). These densities can be expressed in terms of Green functions ^ss'(r, r') 
(s, s' = il). The one-body density of particles of sign s is ns(r) = mgss{'^^ r), the two- 
body truncated densities are rifg^{r,v') = — m^^rss'l^, r')5rg/s(r', r), etc... The matrix g is 
defined by ^ = m~^K{l + K)~^, and therefore g++ and g-+ satisfy the integral equations 

tTTi r \ 1 

5'++(ri,r2) + — / (fr{—— -£f++(r, ra) + — — -6/_+(r,r2)) 

4:W Jd smh/c(2i — z) smhA;(2;i — z) 

i 

4W siiih.k{zi — Z2) 

tTTl f \ 1 

^_+(ri,r2) - — / d^r{ -g++{r,r2) + . -g_+{r,r2)) 

4:W Jd smhA;(2;i — z) smhA;(2;i — z) 

—i 

4:W siiih k{zi — Z2) 

(similar equations hold for g and g-\ ). Differentiating these equations with respect to 

zi and zi respectively, we find differential equations which can be combined into^ 

(m^ - Ai) g++{ri, ra) = m5{ri - ra) 

g-+{ri, ra) = — -— + i— g++{ri, ra) 
m \oxi oyij 

^ Refs. 10 and 11 use slightly different representations of the operator K. In the present 
section, following ref.ll rather than ref. 10, we have in cq. (3.10) a spinor ^ = {i^^ix) 

rather than {i/j, x) and now a, g |_ which is i times the g \. of ref. 10. This has no incidence 

on the physical quantities. 
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Again the boundary conditions are given by the integral equations : 

g++{xi, 0; ra) = -g-+{xi, 0; ra) and g++{xi, W; ra) = g-+ixi, W; ra) 
For the present geometry it is useful to work with the Fourier transforms defined by 

/ + 00 J] 
-oo ^Tl" 

The solution for ^_|_+ is: 

[— ( + Z)e«(2'i+^2)-Kiv _ _ i^^-^(y,+y,)+nw\ _ !!!(^ _ TO)e-«^cosh«(yi - ya)] 
2k, V / K 

X [(« - m)e-'^^ + (av + m)e"^] ~' 



with K = \/P + m^. As ^ +oo for fixed values of j/i and y2, one finds after some 
manipulations 

m f K — I — m \ - w 

g++{yi,y2,l) = — exp[-/t|yi - y2\] + — exp[-K(yi + 7/2)] I + 0(e ) 

Up to an exponentially small term, this expression is the same one as in the case of one ideal 
conductor wall (see eq. (3.24) of ref.lO). There are no algebraic correction in 1/W. This 
confirms our basic assumption of Section 3.2 that the correlation functions are unchanged, 
up to exponentially small effects, by the presence of the second wall. 

3.4. The two-dimensional one-component plasma 

This is a Coulomb system with only one species of mobile particles, of charge q, 
embedded in a background carrying a uniform charge density of the opposite sign. In two 
dimensions, the model is exactly solvable^-*^^'-^^) when = 2. The half-plane geometry, 
with only one ideal conductor wall, has already been studied^^^^. Since the particles repel 
each other, there is no pair collapse and one can deal with point-particles. In presence 
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of ideal conductor walls however, one should prevent the collapse of a particle onto its 
electric image ; here we shall assume that an infinite potential barrier keeps the particles 
away from the walls by a distance £, while the background extends up to the walls (this is 
slightly different from what has been assumed in ref. 15). 

For the present strip geometry, the simplest approach is to take an appropriate limit 
of the two-component case. Instead of one fugacity, one can introduce different fugacities 
for positive and negative particles, and make the negative particle fugacity go to zero. 
Furthermore, a uniform background of charge density —qr] generates an electric potential 
—7rqr]y{W — y) (it has the correct Laplacian 2TTqr} and it vanishes on the walls at y = and 
y = W) ; this potential can be taken into account by replacing the constant fugacity C, of 
the positive particles by a position-dependent fugacity C exp[27r?7y(VF — y)] (here fiq^ = 2), 
and adding to the grand potential the background self-energy which is (l/&)(5~^'nrfW^ 
per unit length. Then, the set of equations (3.8) reduces to one equation 



where D now is the strip between y = s and y = W — e. 

Since z — z' does not vanish if r, r' G -D, (3.17) indicates that ip{r) is an analytical 
function of z. Because of the translational invariance along the a;-axis, ip must depend 
on X through a factor exp(zi;a;), and therefore ip = exp{ivz). Using this ijj in (3.17) and 
performing the integral upon x' gives the relation between v and A 



Using this A in (3.9) (with the background self-energy added), performing the integral 
upon y', and making the change of variable v = 2TT'r]Wt, gives for the grand potential 




(3.17) 




w-e 





6 




m 




+ 2r]W / dt In 



1 + 



2^2^ 2 cosh{27vriWH) 




2 W 




17 



where erf is the error function. Now we spht the integral over t into an integral going from 
to 1/2 and another one going from 1/2 to +oo. We write 



6 



= 2r]W I lne-"^^'(i-2*)'/2dt 
Jo 



and add this term to the first integral. The resulting integral from to 1/2 can be separated 
in three terms which give, in the limit W -\-oo, up to exponentially small terms, the bulk 
pressure, a first contribution for the surface tension and the universal finite-size correction 



(3pi,W = 2Wr] / dt In = rjW In 
Jo 

.1/2 



-2/?7^i) = 2Wr] f dtln 
Jo 



(3.18) 



+ 



m 



erf {V^W{^ _ t __£,)) + erf {V^Wi^ + t - ^)) 



dt 




TT ^0 



+ -(1 + erf {t - ^/2^e))]dt + 0(6"^^) (3.19) 
m 2 



dt ln[e 

iW 'Jo ^ ) \ ) 



2AW 

The integral from 1/2 to +oo contributes only to the surface tension; it can be reexpressed 
as 



-2/37?) = ,/^ r^^dtln 




TT Jo 



1 + ^e* (1-erf {t + e^)) 



+ 0(6-^^) (3.20) 



The final result for the grand potential is 



/3a; = -ppbW + 2/3(7(1) + 7f )) + ^ + 0(e— ^) 



3.5. Strip with plain hard walls 

The solvable two-dimensional models can also be used for studying the case of a strip 
with plain hard walls. The walls confine the particles, but there is no Dirichlet boundary 
condition for the electric potential, which freely "leaks out" . The Coulomb interaction is 
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just the free space one — ln(|r — r'|/a). In that case, there is no finite-size correction of 
order 1/W in the grand potential or free energy per unit length of the strip, as shown 
hereafter. 

For the two-component plasma, the grand partition function is still given by (3.9) and 
(3.11), but the boundary conditions now are^^^ 

- on the boundary y = 0, i/j = go, di/j/dz = ho 

- on the boundary y = W, ip = gw, d'^/dz = hw 

where go, ho (resp. gw, hw) are analytical functions in the domain y < (resp. y > W) 
vanishing at infinity. For i/j of the form (3.12), these boundary conditions become 

^/;{x,0) = , ^{x,W) = {v > 0) 
oz 

dib 

ijj{x, W) ^0 , —^{x,0) = {v <0) . 

dz 

For each value of v, these equations give two linear homogeneous equations for the coeffi- 
cients A and B, the compatibility of which leads to the relation 



cosh{W\/v^ + A2) + |^;| 



smh{WVv^ + A2) 
Vv"^ + A2 



= 



instead of (3.13). Now 



m = e 



-w\ 



1 ^TTT o\ I , sinh(VF\/t'^ + 

cosh(iyv^^2 _^ ^2) _^ 



_j_ 



and, in the limit W ^ oo. 



ln/(-m) ~ W [^Jv'^ + m2) - \v\^ + ^ ( 



+ O (e-™^) 



which gives for the grand potential uj per unit length 



(3uj = 



-PpbW + 2/57 + (e-"^^) 
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The bulk pressure pb is the same as in (3.15), as expected. The surface tension 7 is now 
given by 

(an aheady known result'-^'^''). And there is no finite-size correction, except for an expo- 
nentiaUy smaU one. 

For the one-component plasma in a strip with plain hard walls, the canonical free 
energy has been computed in ref. 16. In our notation, with the de Broglie wavelength and 
the length a in the logarithmic potential taken as unity, the result for the free energy / 
per unit length is 

2 27r2 V27ryo 2 
where Y — W{7r'r]/2)^^'^ and erf is the error function. In the large-M^ limit one finds 

T] rj rn~ f°° l + erf(i) g"'^^^' 

2 27r2 V 27r 7o 2 Att^tj^^W^ 

The first term corresponds to the bulk free energy, the second one to the surface tensions^^^^ 

on both boundaries, and the next term is indeed more than exponentially small ; there is 

no algebraic finite-size correction. 

4. DISK 

This Section is about a two-dimensional Coulomb system in a disk of radius 

R, with an ideal conductor boundary. We derive the expansion (1.4) in general, and check 
it on two solvable models. 

4.1. General derivation 

We have not been able to make a direct derivation of (1.4) using, for instance, screen- 
ing properties of the Coulomb system. We can only rely on the heuristic method of Section 
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2, leading to the relation (2.5) between the Coulomb system and the Gaussian field the- 
ory. For the latter, it has been shown*^^'^'^®^ that — In Zq has a universal finite-size term 
-(1/6) In i?. Prom (2.5), we obtain (1.4). 

4.2. The two-component plasma 

The two-component plasma model, at F = 2, is also solvable in the present geometry. 
In terms of complex coordinates z = re*^, with the origin at the center O of the disk, the 
Coulomb potential solution of (2.3) with Dirichlet boundary conditions is obtained by the 
method of images as 

^(r, r') = - In 



R{z - z') 



ZZ> - i?2 

Like in Section 3.3, we obtain (3.9) for the grand partition function. Now the sum is on 
the eigenvalues A defined by the two coupled integral equations 

27r 
A 

2^ 

where the domain of integration D is the disk. Using the equality 



— =7r5(r-r) 

az z — z' 



we again find the Laplacian equation (3.11) with now the boundary condition x = ^ Wi 
i.e. 

1 .„ / «9 i a\ 

= (4.2) 



r=R 



(again, —A is non-hermitian with that boundary condition). 

Because of the circular symmetry, we look for eigenfunctions of the form 
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where Ii is a modified Bessel function. The boundary condition (4.2) imposes 

I[{\R) - ^leiXR) - h{\R) = 

i.e. 

h+i{\R) - h{\R) = Q . (4.3) 

All the eigenvalues A are obtained by taking, for each value of £ G Z, all the non-zero roots 
of (4.3). Since (4.3) is invariant under the transformation £ —£ — 1, it is enough to 
consider the eigenvalues Xg associated to £ G -CV, counting each of them twice. Then, as in 
Section 3.3, for each value of £ a partial summation of (3.9) can be made by noting that 
the corr responding eigenvalues Xi are the zeros of the function f£{z) : 

/£(^)=^!(^) MzR) - Ie+,{zR)] , e>0 (4.4) 

The prefactor in (4.4) has been chosen in such a way that /g(0) = 1. Therefore, the entire 
function fi{z) can be written as the product 



and 



Eln(l + g^.nn(l + 3=lnM-".) 



making (3.9) a sum on £ only 



lnS = 2^1n|£! (^-^^ [/^(mi?) + /^+i(mi?)] | . (4.5) 

For obtaining a large- asymptotic expansion of (4.5), one can use the Debye expansion 
of IiimR). This expansion is the one which is appropriate here, since both mR and £ be- 
come large. It is convenient to split (4.5) as 



lnS = 51 + ^2 
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where 



and 



£=0 



mR 



l£{mR) 



£=0 



1 + 



l£+l 



(mR) 



himR) 



since is In S for a plain hard wall boundary and has already been evaluated in ref. 2. 
As to S21 rewritten and cut off as 

L 



^2 = 2 J] In 



1 + 



IiimR) mR 



it can be evaluated by similar techniques, using the Debye expansion of In and and using 
for the sum on I the Euler-MacLaurin summation formula . If we now take into account 
the particle-particle and particle-wall short range repulsions by choosing the cutoff as 
L = R/a, we find 



2 1 TT /I 

mi?ln + -ln2 + { 

ma 2 4 \mR 



(4.6) 



In (4.6), all terms which vanish as the cutoff 00 have already been removed. 

Adding Si (from ref. 2) and S2, one finds for the grand potential Q, in the case of an 
ideal conductor boundary. 



pn = -InE 



-(3pbTTR^ + /?7c27ri? + - In mi? - — - - In 2 - 2C'(-1) + - (4.7) 

6 12 6 ^ ^ ^ 4 ^ ^ 



where the bulk pressure pi, and the surface tension 7c are, as expected, given by the same 
expressions (3.15) and (3.16) as in the strip geometry. And (4.7) has the expected universal 
finite-size correction {l/Q)£nR. 

One should note that there is no In R term in S2 ; it comes entirely from S\ . In other 
words, the universal finite-size correction is the same one for a plain hard wall and for an 
ideal conductor wall, unlike in the strip geometry. 
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4.3. The one-component plasma 

We consider the same model as in Section 3.4, but in the disk geometry. The electric 
potential generated by the background is —^TTr]q{R^ — r^). The background self-energy 
adds to the grand potential a term ^P~^(TTr]R'^)'^ The eigenvalue integral equation now is 



XR f exp[7iri{R^ - H)] ^ , 



2n 



i?2 



iP(r')d^r' = V(r) 



(4.8) 



Id z'z 

Again we notice that t/^ is analytical. Because of the rotational invariance, must depend 
on 9 through a factor exp(i£9), and therefore ^^{r) = with £ a non-negative integer. 
Using this form in (4.8) we have the relation between A and £ 

/ ^ \ 2^+1 



A exp[7r?7-R 



— I exp[— 7r?7r^](ir = 1 
RJ 



The grand potential becomes, after some manipulations. 



pn = -{ttvRY 



-l-oo 



e=o 



1 + , 7(^ + 1, ^r]{R - ef) 



(4.9) 



2(7r?7i?2)^+i 

where 7(n, x) = /J^ t"^"^ exp(— i)cii is the incomplete gamma function. Let N = tttjR'^. 
We split the sum over £ in eq. (4.9) into a sum Si for < £ < N and a sum ^'2 for £ > N. 
Using 

f2^ J£-m'^ £-N. 

(4.10) 



gJV^-(^+i)^, ^ J_exp \- + 

V TV ^ 27V 27V ^ 



we have in the limit i? ^ -|-oo 



, mR , 

5i = (iV + l)ln(^) 



JV 



+ ^ [ - (£ + l)lniV + Ar + ln£!] 



N 



+ E 1^ [^e-^ + -Mi + 1, TvviR - ef)] 



(4.11) 
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Using Stirling's formula , the asymptotic formulas for i? ^ +00 

^ln£! = (Ar + l)lnAr! IniV 



iV2 1 1 

+ N\nN In AT + 0(1) 

4 2 12 ^ ^ 



(4.12) 



and 



we obtain, in the small -e limit, 



N -e , — 

1 + erf ( ^ - e^/27^r] 



^2N 



+ o{ 



(4.13) 



1 / m \ 1 



+ J2^R I dt\n T-^e"* + -(1 + erf {t - eV2Tr eta))] + 0(1) (4.14) 
Jo m 2 

In ^2 we use (4.13) and (4.10) to find 



S2 = ^/2^R / dt In 
Jo 



m ^2 

1 + 



2v^ 



erf {t + ey/2^)^ 



(4.15) 



Putting (4.14) and (4.15) in (4.9) we have 



= -PpbTTR' + 27ri?/3(7^^) + 7^2)) + i In 



m 



{7rr])'/^R 



+ 0(1) (4.16) 



where the bulk pressure pi, is given by (3.18) and 7^^^ and 7^^^ are given by (3.19) and 
(3.20) in the small-e limit. We find the universal finite-size correction | Ini?. 

It is interesting to compare the surface tension for the one-component plasma in the 
limit e — > and the surface tension for the two-component plasma in the limit a — > 0. 
We have in that limit for the one-component plasma 



m 



47r 



7c"' ~ — In e 



and for the two-component plasma 



m 

~ — In a 
27r 
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The surface tension of the one-component plasma is one-half the one of the two-component 
plasma. 

5. CONCLUSION 

The grand potential of a Coulomb system has universal finite-size corrections very- 
similar to the ones which occur for a simple critical system : the massless Gaussian field 
theory. The boundary conditions however play an important role. In the geometry of a 
slab, if the boundaries are ideal conductor walls, a Coulomb system has a grand potential 
with a universal finite-size correction which is just the opposite of the one found for a 
Gaussian field with Dirichlet boundary conditions; but if the boundaries are plain hard 
walls, the Coulomb system no longer exhibits that finite-size correction. In the geometry of 
a disk, a two-dimensional Coulomb system has a finite-size correction to its grand potential 
which is the same one for an ideal conductor wall and for the previously studied*^^^ case of 
a plain hard wall (in both cases, again the correction is opposite to the one occuring for the 
Gaussian field). An argument given in ref. 2 essentially says that a Coulomb system with 
plain hard walls looks somewhat like an ideal conductor when seen from the outside; this 
might be relevant for explaining that the curved boundary of a disk generates the same 
correction (1/6) InR in both plain hard wall and ideal conductor wall cases. Nevertheless, 
a unified treatment of the different possible boundary conditions still has to be found. 

In the case of a slab with ideal conductor walls, the finite-size correction to the grand 
potential of a Coulomb system was derived on the sole basis of the screening effect (Section 
3.2). A more general derivation on the basis of the screening effect should be found for 
other geometries and boundary conditions. 
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